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ABSTRACT 


SRI  International  has  completed  a  program  of  research  on  an  innovative  time-differential 
phase  sensor  suitable  for  detecting  the  fast  changing  structure  within  airborne  turbulence.  This 
innovation  is  based  on  an  original  incorporation  of  time-domain-stimulated  photon  echo  into  a 
phase-conjugate-mirror  of  an  optical  novelty  filter  setup.  The  unique  advantage  of  this  photon 
echo  novelty  filter  is  that  its  fast  response  time  is  completely  programmable  through  external 
excitation  timing.  Therefore,  it  has  an  outstanding  potential  for  tracking  fast  changing  phase 
distortions  in  real  time. 

Using  Eu  and  Pr  doped  yttrium  silicate  crystal  as  the  hole-burning  materials  for  the 
stimulated  photon  echo,  we  demonstrated  the  concept  and  basic  performance  of  the  photon  echo 
novelty  filter.  The  sensitivity  of  the  filter  to  a  sudden  change  of  optical  path  length  is  measured 
to  be  better  than  0.05  wavelength  (~0.6  jjm).  The  fast  temporal  tracking  ability  of  the  filter  is 
proved  by  its  accurate  measurement  of  phase  modulation  events  from  a  nanosecond  through 
microsecond  time  scale.  The  filter's  holographic  memory  and  time-differential  processing  of 
nonuniform  phase-only  distortion  is  proved.  Our  results  show  that  the  photon  echo  novelty  filter 
is  well  worth  further  investigation  and  development  into  a  powerful  sensor  of  the  temporal  and 
spatial  stmctures  of  airborne  turbulence. 
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INTRODUCTION 


Airborne  turbulence  has  been  a  long-standing  problem  in  aero-optics  [1].  Optical 
imaging  and  the  laser  projection  wavefront  from  a  flying  airplane  platform  are  strongly  distorted 
by  the  surrounding  turbulent  air  layer  that  the  light  traverses.  This  problem  becomes  much  more 
severe  for  rear  looking  directions  [2]  and  is  aggravated  by  the  turbulence’s  fast  change  and  short 
coherence  time  on  microsecond  time  scale  [3].  The  coherence  time  is  three  orders  of  magnitude 
shorter  than  the  response  time  of  current  adaptive  optics  [4].  A  need  for  a  correspondingly  fast 
phase  sensor  and  compensator  is  thus  evident.  Because  little  is  known  about  the  temporal  and 
spatial  structures  inside  airborne  turbulence,  a  reasonable  approach  to  the  problem  is  to  study  the 
turbulence  using  a  fast  time-differential  phase  sensor  that  can  temporally  resolve  the  fast  change. 

An  optical  novelty  filter  is  a  suitable  tool  for  resolving  the  fast-changing  feature  in  an 
optical  input  [5].  Conventional  optical  novelty  filters  are  based  mainly  on  photorefractive  crystals, 
which  take  milliseconds  to  react  to  a  constantly  changing  input  [6].  Clearly,  a  new  novelty  filter  is 
needed  to  react  to  the  fast  and  constantly  changing  input  through  airborne  turbulence. 

This  report  describes  a  new  optical  novelty  filter  based  on  spectral  hole-burning  crystals 
through  time-domain  stimulated  photon  echo.  This  new  filter  is  unique  in  its  programmable 
response  time,  which  can  be  changed  from  nanoseconds  to  microseconds. 

Stimulated  photon  echo  technologies  have  been  developed  by  SRI  International  for  many 
years  [7-9].  Our  good  understanding  of  the  basic  physics  and  optical  memory  application  of 
stimulated  photon  echo  in  spectral  hole-burning  crystals  encouraged  us  to  construct  a  phase-only 
optical  memory  and  a  phase-conjugate-mirror  installed  in  the  new  novelty  filter. 

During  the  past  two  years,  we  demonstrated  the  concept  and  basic  performances  of  the 
photon  echo  novelty  filter  to  detect  fast  phase  changes.  We  showed  that  the  filter  is  fully  capable 
of  real  time  resolving  a  fast  phase  change  from  a  nanosecond  to  microsecond  time  scale  [10,1 1]. 
The  filter  was  found  to  be  more  suitable  for  detecting  a  fast  phase  change  within  20  (is  than  for 
slower  ones  because  of  its  intrinsic  high-pass  feature.  We  suppressed  the  background  of  the  filter 
to  an  equivalent  phase  change  below  0.05  wavelength  (~0.6  (im).  The  accuracy  and  resolution  of 
the  filter  were  found  to  be  0.02  wavelength.  The  filter’s  holographic  phase-only  memory  is 
proved  by  the  phase -conjugate  recalling  of  a  random  phase  map  [11].  These  results  provide  the 
basis  for  future  development  of  the  filter  into  a  real  time  tracking  sensor  of  fast  phase  change 
through  airborne  turbulence. 
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BACKGROUND 


AIRBORNE  TURBULENCE 

For  an  airplane  flying  at  the  speed  of  sound,  the  air  surrounding  the  plane  is  dragged  by 
the  plane  to  move  forward  at  a  similar  speed.  However,  a  boundaiy  layer  of  air  adheres  to  the 
front  of  the  plane's  body.  The  flow  inside  this  boundary  layer  is  nearly  laminar  and  does  not 
distort  a  light  wave  going  through  it  very  much  [2],  partly  because  the  air  in  the  boundary  layer 
has  a  very  small  speed  transverse  to  a  forward  going  light  wave.  The  wind  speed  longitudinal  to 
the  optical  wave  produces  no  wavefront  distortion.  Therefore,  the  Strehl  ratio  is  quite  high  in  the 
forward  direction  (azimuthal  angle  close  to  0°). 

For  a  light  wave  going  at  an  azimuthal  angle  substantially  larger  than  zero  but  less  than 
90°,  the  air  in  the  boundary  layer  has  a  large  transverse  speed  across  its  wavefront. 

Consequently,  the  wavefront  is  distorted  and  the  Strehl  ratio  is  smaller.  However,  because  the 
flow  in  the  boundary  layer  is  laminar,  the  air  density  and  the  index  of  refraction  do  not  fluctuate 
very  much.  Therefore,  the  wavefront  is  not  severely  distorted  at  the  front  of  the  plane  and  the 
Strehl  ratio  is  tolerable. 

The  picture  changes  completely  at  the  back  of  the  airplane  (see  Figure  1).  The  boundary 
layer  begins  to  be  sepai'ated  from  the  body  of  the  plane  at  an  azimuthal  angle  of  90°,  and  this 
separation  is  more  prominent  at  larger  azimuthal  angles  [2].  In  the  latter  region,  the  boundary 
layer  becomes  a  completely  separated  flow.  This  flow  still  moving  at  a  high  speed  with  the 
plane,  is  surrounded  by  the  ambient  air,  which  is  only  drifting  slowly.  The  large  velocity 
difference  at  the  interface  between  the  separated  flow  and  the  ambient  air  creates  a  shear  layer 
that  is  fully  turbulent  [2],  and  the  strength  of  turbulence  is  very  strong  in  that  region.  The  strong 
fluctuations  of  air  density  cause  strong  fluctuations  of  the  index  of  refraction. 

To  make  things  worse,  the  shear  layer  is  projected  over  a  long  distance  behind  the  plane. 
An  optical  beam  propagating  at  large  azimuthal  angle  inevitably  travels  a  long  path  length 
through  this  shear  layer,  which  accumulates  into  great  wavefront  distortions.  For  the  same  path 
length,  the  corresponding  wavefront  distortion  is  less  serious  for  longer  optical  wavelength. 
However,  shorter  wavelength  is  favored  in  applications  requiring  tighter  diffraction  limited 
focus.  To  show  how  severe  the  problem  is,  it  has  been  reported  that  the  phase  distortion  at  the 
back  of  an  airplane  turret  can  approach  nearly  two  wavelengths  (-0.7  pm)  [2]. 
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Figure  1. 


The  other  harsh  character  of  airborne  turbulence  is  its  fast  change  and  short  coherence 
time.  The  phase  distortion  is  reported  to  rise  to  80%  of  its  full-scale  value  within  approximately 
200  |J,s.  In  other  words,  the  turbulent  phase  change  is  coherently  trackable  only  within  about 
200  )is.  To  correct  this  phase  distortion,  an  adaptive  system  needs  to  measure  the  phase  and 
initiate  a  compensation  within  that  coherence  time.  This  time  is  much  shorter  than  the  response 
time  of  Shack-Hartmann  sensor  and  deformable  mirrors  within  the  current  adaptive  optics 
technologies  [4]. 

Such  high  bandwidth  is  also  predicted  in  theory  by  the  Greenwood  frequency  [12],  which 
is  defined  as  the  bandwidth  required  to  fully  compensate  for  wavefront  distortion.  When  the 
transverse  wind  velocity  is  constant  along  the  optical  path,  the  Greenwood  frequency  can  be 
written  as 


fc  =  0.426  v/ro 

where  v  is  the  wind  velocity  and  ro  is  the  Fried  parameter  [12]  for  the  turbulence;  v  is 
approximately  200  m/s  for  the  wind  around  an  airplane  flying  at  a  Mach  number  of  one.  ro  is 
defined  as 


ro  =  0.185(X^/Cn^Az)^^^ 
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where  X  is  the  optical  wavelength,  Cn^  is  the  refractive  index  structure  constant,  which  tells  the 
strength  of  fluctuations  of  the  index  of  refraction,  and  Az  is  the  length  of  the  optical  path,  or  the 
thickness  of  the  turbulence  layer. 

Cn^  is  defined  for  an  isotropic  turbulence  predictable  by  the  Kolmogorov  theory: 

<[n(ri  +  r)  -  n(ri)]^>  =  Cn^ 


where  n  is  the  index  of  refraction,  <  >  means  to  average  over  ri,  and  r  =  Irl. 

According  to  a  wind  tunnel  experiment  [2],  the  average  phase  variation  across  an  optical 
aperture  of  10  cm  is  about  one  wave  (-0.5  |a.m).  This  phase  variation  can  be  approximated  as 

Af  =  [<[n(ri  +  r)  -  n(ri)]^>]  Az 


Using  the  aperture  r  =  10  cm  and  the  experimental  finding  of  shear  layer  thickness  Az  =  10  cm 
[2],  Cn^  =  5  *  10'^  which  is  significantly  larger  than  that  for  natural  atmosphere 
unperturbed  by  aircraft  [13],  about  10'^^  to  10*^^  From  Cn^,  iD  =  2  mm.  From 
V  =  200  m/s,  fG  =  40  kHz,  which  shows  that  the  bandwidth  of  the  airborne  turbulence  is  on  the 
order  of  100  kHz.  Current  adaptive  optics,  however,  has  a  bandwidth  of  a  few  kilohertz  [4,14]. 

Despite  the  complexity  of  airborne  turbulence,  some  coherent  structures  may  exist  inside 
the  shear  layer  [15].  Because  these  structures  may  be  dominant  in  causing  large-scale  air  density 
fluctuations,  we  urgently  need  to  understand  their  temporal  and  spatial  distributions.  This 
knowledge  will  certainly  provide  a  guide  to  developing  fast  adaptive  sensors  and  compensators, 
the  ultimate  goals. 


OPTICAL  NOVELTY  FILTER 

An  optical  novelty  filter  [5]  is  a  high-pass  filter  that  passes  rapidly  changing  input  light 
and  rejects  slowly  changing  input  light.  Many  reported  novelty  filters  [6,10,16-18]  are  in  fact 
time-differential  filters.  By  optically  subtracting  the  current  input  light  from  the  past,  these 
filters  allow  only  the  new  features  in  the  input  to  go  to  the  output.  The  interval  between  the 
current  and  the  past  is  determined  by  the  filter's  response  time. 


4 


A  novelty  filter  has  a  finite  response  time  so  that  it  responds  to  the  past  input  at  the 
current  time.  The  filter  is  often  set  up  as  a  self-aligning  interferometer  [5].  It  consists  of  a 
beamsplitter  at  the  entrance  and  a  phase-conjugate  mirror  at  the  end  (see  Figure  2).  The 
beamsplitter  divides  the  input  light  into  two  copies  along  two  arms,  both  of  them  incident  on  the 
conjugate  mirror. 


o(t) 


The  phase-conjugate  mirror  has  a  finite  response  time.  Instead  of  reflecting  the  incident 
light  instantaneously,  the  mirror  first  stores  the  incident  light  amplitude  and  phase.  After  a 
certain  response  time,  the  mirror  recalls  the  memory  and  retroreflects  the  incident  light  with  a 
conjugate  phase.  The  reflected  light  then  retraces  the  path  of  the  incident  light  back  to  the 
beamsplitter.  Because  of  the  conjugate  phase,  the  reflected  light  cancels  all  the  past  distortions 
caused  by  the  medium  and  picks  up  the  new  phase.  Therefore,  when  it  reaches  the  end,  the 
reflected  light  contains  a  time -differential  phase  distortion  of  the  medium. 

In  the  novelty  filter,  one  arm  traverses  the  phase  changing  medium  and  the  other  arm 
serves  as  constant  reference.  The  light  waves  along  the  two  arms  recombine  at  the  beamsplitter 
and  destructively  interfere  on  the  exit  side  [5].  The  time-differential  phase  of  the  medium  is  thus 
presented  as  the  interference  signal  at  the  output. 

Phase-conjugation  has  been  observed  by  either  stimulated  Brillouin  scattering  or  four- 
wave-mixing  [19].  Some  well-known  materials  for  four-wave-mixing  are  photorefractive 
crystals  and  Kerr  media  [20].  Although  these  materials  are  highly  efficient  and  easy  to  use,  their 
response  times  are  tightly  constrained  by  material  characteristics,  such  as  diffusion  time  [21]  or 
relaxation  time  [22].  Therefore,  they  cannot  provide  a  flexible  response  time  in  accord  with  the 
coherence  time  of  turbulence  structures. 
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TIME-DOMAIN  STIMULATED  PHOTON  ECHO 


In  a  three-pulse  stimulated  photon  echo  experiment  [23],  the  first  two  laser  pulses  burn  a 
persistent  spectral  hole  pattern  in  the  ground  state  population  of  a  material,  such  as  a  rare  earth 
doped  crystal  at  low  temperature.  After  a  long  waiting  time,  the  third  laser  pulse  excites  the 
material  again  and  stimulates  a  coherent  radiation  pulse  that  is  delayed  in  time.  This  time  delay 
is  equal  to  the  time  delay  between  the  first  two  laser  pulses,  hence  the  name  echo.  In  other 
words,  the  initial  time  delay  between  the  first  two  laser  pulses,  is  memorized  by  the  material  and 
recalled  exactly  through  the  echo’s  timing.  This  provides  the  basis  of  optical  memory  of 
temporal  modulation. 

However,  the  spectral  hole  pattern  burnt  in  the  material  by  the  first  two  laser  pulses  is  a 
spectral  hologram.  Although  not  overlapping  in  time,  the  two  laser  pulses  overlap  in  frequency. 
One  laser  serves  as  the  object  light  and  the  other  as  the  reference  [24].  The  third  laser  pulse 
serves  as  a  reference  again  to  read  out  the  echo  as  a  holographic  image  of  the  object.  This 
provides  the  basis  of  optical  memory  of  spatial  modulation.  Therefore,  an  optical  memory  can 
be  installed  either  through  temporal  modulation  of  the  first  or  second  laser  pulse  shape,  or 
through  spatial  modulation  of  these  pulses’  intensity  or  phase  profile  [8,9]. 

Spectral  hologram  needs  to  be  explained.  Very  narrow  spectral  holes  can  be  burnt  in  the 
absorption  spectrum  of  a  rare  earth  doped  crystal  at  low  temperature.  This  indicates  that  each 
rare  earth  ion  has  a  long  coherence  time.  The  spectral  holes  burnt  by  the  combined  first  two  laser 
pulses  are  caused  by  the  interference  between  them.  Although  the  two  laser  pulses  do  not 
overlap  in  time,  the  rare  earth  ions  excited  by  the  first  pulse  oscillate  coherently  after  that  pulse  is 
gone  and  react  to  the  coming  second  pulse.  Therefore,  these  ions  provide  the  basis  of  coherent 
interference  between  the  two  laser  pulses.  If  each  laser  pulse  is  weak  in  intensity,  the  result  of 
this  interference  is  a  second-order  spectral  modulation  of  the  ions'  ground  state  population.  This 
modulation  is  a  long-lived  spectral  hole  pattern  or  a  spectral  hologram. 

After  a  long  waiting  time,  the  third  laser  pulse  excites  the  rare  earth  ions  from  this 
modulated  ground  state  and  stimulates  their  coherent  dipole  radiations.  These  ions'  radiations 
become  aligned  in  phase  at  a  delayed  time,  which  is  equal  to  the  time  delay  between  the  first  two 
laser  pulses  [23].  Thus  the  macroscopic  radiation  reaches  maximum  at  that  time  in  the 
appearance  of  an  echo  pulse. 
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In  general,  the  echo  in  the  frequency  space  can  be  viewed  as  the  linear  diffraction  of  the 
third  laser  pulse  from  the  spectral  hologram  formed  by  the  coherent  interference  between  the  first 
and  second  laser  pulses  in  a  material  with  narrow  homogeneous  linewidth: 


Ee(r,a))  ~  E3(r,co){E2(r,a))[  Ei(r,co)]*}  (1) 

where  El ,2,3  (co)  are  the  amplitudes  of  the  three  laser  pulses  in  the  frequency  space  [25],  and  the 
bracket  {  }  contains  the  spectral  hologram’s  amplitude.  Eq.  (1)  is  valid  for  weak  laser  intensities. 

It  is  then  evident  that  the  echo  is  a  conjugate  image  of  Ei(r,(o)  when  E3(r,co)  E2(r,(0)  is 
uniform  in  space,  or  a  real  image  of  E2(r,(0)  if  E3(r,0})[  Ei(r,(jD)]*  is  uniform  in  space. 

PHASE-CONJUGATE  REFLECTION 

When  E3(r,co)  E2(r,Q))  is  uniform  in  space,  Ee(r,C0)~[  Ei(r,co)]*.  The  radiation  from  a 
material  is  the  complex  conjugate  of  an  incident  laser.  This  is  the  first  feature  of  phase  conjugate 
reflection.  The  echo  is  pointed  at  ke  =  k3  +  k2  -  ki,  where  ki,  k2,  and  k3  are  the  wave  vectors  of 
the  three  laser  pulses  in  the  nonlinear  material.  A  backward  photon  echo  is  established  when  the 
second  and  third  laser  pulses  are  antiparallel,  i.e.,  k3  -f  k2  =  0.  Under  this  condition,  ke  =  -ki. 
The  new  radiation  pointing  as  a  retroflection  of  the  incident  laser  is  the  second  feature  of  phase 
conjugate  reflection.  Therefore,  the  backward  photon  echo  is  a  phase-conjugate  reflection  of  the 
first  laser  pulse.  The  backwai'd  configuration  is  similar  to  most  four-wave-mixing  phase- 
conjugate -mirrors  [22]. 

However,  the  photon  echo  phase-conjugate-mirror  (PCM)  has  a  unique  character  that 
other  PCMs  do  not  have:  the  arbitrary  time  delay  of  its  phase-conjugate  reflection.  The  echo  is 
radiated  at  t  =  t3  t2  -  ti,  where  ti,  t2,  and  t3  are  the  incidence  times  of  the  three  laser  pulses  in 
the  nonlinear  material.  From  the  incidence  time  of  the  first  laser  pulse,  the  echo  is  delayed  until 
(t2  -  ti)  after  the  third  laser  pulse  anives.  Therefore,  one  can  delay  the  phase-conjugate  reflection 
arbitrarily  by  delaying  the  third  laser  pulse,  as  long  as  the  spectral  hologram  is  long-lived. 
Because  of  the  ease  of  controlling  the  echo's  timing  by  programming  external  excitation  delay, 
the  photon  echo  PCM  has  a  very  flexible  response  time.  This  conveniently  suits  the  study  of 
dynamic  phase  change  on  any  fast  time  scale  from  nanosecond  to  microsecond. 

The  reflectivity  of  the  photon  echo  PCM  depends  on  several  parameters.  The  first  laser 
pulse  should  be  aligned  to  strike  the  PCM  on  the  effective  area  that  is  pumped  by  the  second  and 
the  third  laser  pulses.  If  the  pump  laser  pulses  have  a  Gaussian  spatial  profile,  then  clearly  the 
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reflectivity  at  the  center  of  the  pumped  area  is  different  from  the  wing.  The  reflectivity  also 
depends  on  the  time  delay  between  the  laser  pulses.  The  delay  between  the  first  and  second 
pulses  should  be  shorter  than  the  coherence  time  of  the  rare  earth  ions  excited  in  the  PCM. 
Otherwise,  there  is  no  reflectivity  because  no  spectral  hologram  can  be  burnt.  This  requirement 
is  similar  to  the  coherence  requirement  in  conventional  holography  [26,27]. 

OUTPUT  OF  PHOTON  ECHO  NOVELTY  FILTER 

In  the  photon  echo  novelty  filter,  the  first  laser  pulse  is  split  by  a  beamsplitter  into  two 
beams  along  two  paths.  The  two  beams  are  finally  incident  on  the  rare  earth  doped  crystal  at  two 
ki  directions  (see  Figure  3).  Because  they  come  along  two  separate  paths,  these  two  beams 
differ  in  their  incidence  direction  on  the  phase-conjugate-mirror.  If  these  two  beams  strike  the 
same  effective  area  of  the  phase-conjugate-mirror,  the  reflectivities  for  them  are  the  same 
(assuming  they  are  coincident).  Therefore,  when  either  arm  of  the  novelty  filter  is  blocked,  we 
should  measure  the  same  echo  intensity  along  the  other  arm  at  the  output  of  the  filter  because  the 
two  arms  are  symmetric  about  the  beamsplitter.  That  is,  the  light  on  either  arm  experiences  one 
reflection  and  one  transmission  at  the  beamsplitter  to  go  from  the  input  to  the  output  of  the  filter. 

If  the  echo  along  one  arm  sees  a  new  phase  change  of  5  while  the  other  arm  remains 
constant,  the  output  of  the  novelty  filter  (channel  y)  is 

I(y)  =  2(l -cos6)I(e)  (2) 

where  1(e)  is  the  output  when  either  arm  is  blocked.  Here  we  assume  that  the  two  ki's  have  no 
cross  talk  in  the  nonlinear  crystal  so  that  blocking  one  does  not  affect  the  other.  This  is  true  when 
both  two  ki's  are  weak  in  intensity  so  that  their  combination  in  the  crystal  is  also  weak  in  intensity. 

Using  Eq.  (2),  we  can  find  6  by  measuring  I(y)  and  1(e).  These  measurements  are  not 
done  at  the  same  time  because  one  arm  needs  to  be  blocked  to  measure  1(e).  To  find  5  without 
measuring  1(e),  we  refer  to  the  monotony  output  [5]  at  the  beamsplitter.  This  output  is  the 
recombined  echo  that  goes  back  to  the  light  source.  Because  the  monotony  output  is  also 
proportional  to  1(e),  we  can  find  5  by  measuring  the  novelty  (channel  y)  and  the  monotony 
(channel  x)  output  at  the  same  time: 

I(y)/I(x)  =  (1-  cos5  )/(p  +  q  cos5 )  (3) 

where  p  and  q  are  constants  defined  in  appendix  A. 


8 


EXPERIMENTAL  CONDITION 


PROOF  OF  CONCEPT  ON  NANOSECOND  TIME  SCALE 

The  laser  pulses  for  the  experiment  were  produced  by  a  dye  laser  (SpectraPhysics  PDL-1) 
pumped  by  the  second  harmonic  of  a  Q-switched  Nd:YAG  laser  (Quanta  Ray  DCR-1).  Running 
at  10  Hz  repetition  rate,  the  laser  pulse  was  7  ns  in  duration,  6  GHz  in  bandwidth,  and  580  nm  in 
wavelength.  This  wavelength  was  tuned  to  excite  the  ^Fq-^Dq  transition  of  a  Eu-doped  crystal 
that  was  used  for  the  phase-conjugate-mirror.  The  laser  was  linearly  polarized  by  a  Gian  prism. 

The  laser  pulse  was  split  into  three  pulses  at  times  ti,  t2,  and  t3  (see  Figure  4).  A  white 
cell  introduced  delays  of  about  24  ns  and  1 14  ns  for  t2  -  ti  and  t3  -ti,  respectively.  The  first  laser 
pulse  was  split  into  two  copies  that  were  both  focused  to  0.2  mm  on  the  europium  crystal.  They 
were  almost  coincident.  The  second  and  the  third  laser  pulses  were  spatially  filtered  and 
collimated.  They  were  aligned  to  overlap  and  counterpropagate  through  the  europium  crystal 
with  the  same  beam  diameter  of  1  mm.  The  angle  between  either  ki  and  k3  was  about  50  mrad. 
When  the  intensity  of  any  one  of  three  laser  pulses  was  attenuated  by  50%,  the  echo  intensity 
also  dropped  by  50%.  This  indicated  that  the  laser  intensity  is  weak  enough. 

The  europium  doped  crystal  is  0.1%  Eu^+:Y2Si05  (Scientific  Materials),  which  is  8  mm 
in  diameter  and  7  mm  long.  It  was  mounted  in  a  closed-cycle  refrigerator  (RMC  model  22).  The 
temperature  was  kept  at  19-20  K.  This  temperature  was  the  lowest  point  to  remove  the 
accumulated  photon  echo  from  the  crystal  [28].  All  the  echo  data  were  taken  at  site  2  of  the 
crystal  [29].  The  optical  density  of  our  crystal  at  site  2  was  less  than  one. 

An  avalanche  transistor  driver  provided  microsecond-long  high  voltage  pulses  to  a 
Pockels  cell  (KD*P)  installed  along  one  arm  of  the  novelty  filter.  The  KD*P  crystal  generates  a 
phase  modulation  pulse  that  is  linearly  proportional  to  the  high  voltage  pulse.  The  KD*P 
crystal's  response  time  is  a  few  nanoseconds.  The  phase  modulation  pulse  has  a  rise  time  of  a 
few  nanoseconds  and  a  length  of  a  few  microseconds.  This  pulse  was  seen  as  a  step  function  by 
the  nanosecond  light  pulses.  The  beginning  time  of  the  high  voltage  pulse  was  after  the  first 
laser  pulse  but  before  the  echo.  To  use  the  KD*P  as  a  pure  phase  modulator,  we  aligned  one  of 
its  extinction  axes  parallel  to  the  laser's  linear  polarization. 
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laser 


detector 


The  phase  modulations  (6)  along  the  two  extinction  axes  of  KD*P  are  equal  but  opposite 
in  sign  [30].  To  calibrate  5,  we  aligned  the  extinction  axes  at  exactly  45°  from  the  laser's 
polarization.  The  beginning  time  of  the  high  voltage  pulse  was  changed  to  be  before  the  first 
laser  pulse.  We  measured  the  output  intensity  of  the  Pockels  cell  after  a  cross  polarizer.  In  this 
configuration,  Iq  =  Ii(sin5)^,  where  li  and  Iq  were  the  input  and  output  intensity,  respectively. 
Thus  5  was  calibrated.  The  maximum  phase  modulation  in  one  single  KD*P  was  less  than  90°. 
When  we  wanted  a  modulation  more  than  90°,  we  used  two  KD*P  crystals  in  series. 

We  tested  the  filter's  output  under  several  5  values.  We  varied  5  in  two  ways.  First,  we 
varied  it  among  a  set  of  discrete  values  by  connecting  cables  of  different  lengths  to  the  avalanche 
transistor  driver  output.  This  changed  the  capacitive  load  to  the  driver  and  amplitude  of  the  high 
voltage  pulse.  Second,  the  driver  electronics  was  modified  to  shorten  the  high  voltage  pulse  to 
100  ns.  In  this  way,  the  phase  modulation  seen  by  the  echo  was  varied  by  shifting  the  beginning 
time  of  the  high  voltage  pulse  relative  to  the  echo.  Thus  the  echo  saw  a  different  amplitude  of 
phase  change  along  the  100-ns-long  pulse.  To  prevent  the  first  laser  pulse  from  also  seeing  the 
phase  change,  the  high  voltage  pulse  had  to  be  shorter  than  te  -  ti,  where  te  and  ti  were  the  time 
of  the  echo  and  the  first  laser  pulse,  respectively. 

A  light  switch  consisting  of  a  Pockels  cell  and  two  cross  polarizers  was  installed  to 
remove  the  laser  scattering  from  the  echo  at  the  output.  The  echo  was  first  detected  by  a  CCD 
camera.  We  aligned  the  two  ki's  to  achieve  a  Gaussian  spatial  profile  for  the  echoes  and  to 
minimize  the  background  when  the  echoes  canceled  each  other.  The  echoes  at  channel  y  and  x 
were  then  detected  by  fast  photodiodes.  The  average  echo  gave  a  voltage  pulse  of  about  100  mV 
in  the  photodiode.  The  maximum  echo  was  about  10  times  larger  in  intensity.  The  output  of  the 
photodiodes  was  integrated  by  boxcar  integrators  and  then  digitized  by  an  A/D  board  in  a  PC. 

TEST  OF  PERFORMANCE  ON  MICROSECOND  TIME  SCALE 

We  used  two  Q-switched  Nd:YAG  lasers  to  pump  two  pulsed  dye  lasers.  The  dye  lasers 
had  7-ns  pulse  duration,  6-GHz  bandwidth,  and  10-Hz  repetition  rate.  One  dye  laser  provided 
both  the  first  and  second  laser  pulses,  which  were  optically  separated  by  20  ns.  The  other  dye 
laser  provided  the  third  laser  pulse,  which  could  be  continuously  delayed  from  the  first  dye  laser 
(see  Figure  5).  We  chose  either  Eu^’*';Y2Si05  (8-mm  diameter,  7-mm  thickness)  or  Pr^''';Y2Si05 
(8-mm  diameter,  2-mm  thickness)  as  the  nonlinear  material  for  photon  echo.  The  Eu  and  Pr 
crystals  were  held  at  18  and  9  K,  respectively,  in  a  closed  cycle  refrigerator.  At  these 
temperatures,  the  accumulated  photon  echoes  vanished  in  these  crystals. 
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The  dye  lasers  were  tuned  to  excite  either  the  ^Fq-^Do  transition  of  Y2Si05  or  the 
^H4-*D2  transition  of  Pr^'*’:Y2Si05.  Because  there  are  two  sites  in  the  Y2Si05,  the  two  lasers 
were  tuned  to  excite  the  rare  earth  ions  at  the  same  site.  We  used  site  2  for  Eu  and  site  1  for  Pr  to 
take  advantage  of  the  ions'  larger  oscillator  strengths  at  these  sites  [31,32].  The  inhomogeneous 
linewidths  of  Eu  and  Pr  at  these  sites  were  measured  to  be  11  and  49  GHz  (broadened  by  the 
6-GHz  laser  bandwidth),  respectively.  The  peak  absorptions  of  the  laser  in  our  Eu  and  Pr  doped 
crystals  were  18%  and  70%,  respectively. 

Because  Pr  has  about  ten  times  greater  oscillator  strength  than  Eu,  the  echo  signal  from  Pr 
doped  crystal  was  much  larger.  However,  the  laser  intensity  could  also  be  too  high  for  Pr 
because  of  its  large  oscillator  strength.  A  strong  incidence  laser  intensity  could  violate  the 
condition  of  phase-conjugate-reflection.  In  our  experiments,  we  did  not  focus  the  first  laser  pulse 
on  the  Pr  doped  crystal.  In  this  way  we  prevented  the  laser  intensity  from  going  too  high.  For 
both  the  Eu  and  Pr  doped  crystals,  the  optimized  echo  signal  was  bright  enough  to  be  seen  and 
gave  more  than  1  V  of  response  on  fast  photodiodes. 

Two  transparent  media  that  can  generate  phase  changes  were  used  to  test  the  novelty 
filter.  The  first  medium  was  a  KD*P  crystal  to  which  we  supplied  microsecond-long  high 
voltage  pulses.  The  beginning  time  of  the  high  voltage  pulse  was  fixed  relative  to  the  first  two 
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laser  pulses  from  one  dye  laser.  The  third  laser  pulse  (from  the  other  dye  laser)  was  scanned  in 
time  relative  to  the  high  voltage  pulse.  In  this  way,  the  echo  was  scanned  across  the 
microsecond-long  pulse  of  phase  change. 

The  second  medium  was  a  helium  jet.  Compressed  helium  gas  was  ejected  vertically 
from  a  2-mm-diameter  round  nozzle  at  a  speed  about  100  ra/s.  The  horizontal  laser  beam 
intersected  the  jet  above  the  nozzle.  We  selected  helium  because  it  has  a  very  different  index  of 
refraction  from  air  (at  standard  temperature  and  pressure,  n(He)  =  1.000036,  n(air)  =  1.000293). 
Therefore,  a  small  change  in  the  helium  jet  thickness  can  induce  a  large  change  in  the  optical 
path  length  (a  1-mm  change  in  helium  thickness  gives  about  0.25-|im  change  in  optical  path 
length).  We  first  used  fast  photodiodes  to  spatially  integrate  the  echo  intensity  through  the  jet. 
When  the  helium  jet  needed  to  be  spatially  resolved,  it  was  first  imaged  on  the  crystal,  using  the 
first  laser  pulse,  and  then  imaged  on  a  CCD  camera,  using  the  echo. 
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RESULTS  AND  DISCUSSION 


PROOF  OF  CONCEPT  ON  NANOSECOND  TIIVIE  SCALE 
Data  for  Novelty  Filter  Alone 

The  value  of  5  (novelty  filter)  was  calculated  from  I(y)  =  2(1  -  cos§)I(e).  For  each 
measurement  of  1(e)  and  I(y),  we  took  500  laser  shots  and  averaged  the  signal.  1(e)  was 
measured  separately  along  the  two  arms,  and  we  checked  that  they  did  not  differ  more  than  10%. 
The  average  of  1(e)  along  the  two  arms  was  used  in  the  calculation.  We  repeated  this  procedure 
5  times  and  averaged  the  result  to  get  the  final  value  of  5.  a(5)  is  the  standard  deviation  of  5 
from  the  average  in  these  5  measurements. 

5  (calibrated)  was  calculated  from  Iq  =  Ii(sin5)^.  The  error  in  calibration  was  assumed  to 
be  negligible.  The  results  are  given  below. 
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Data  for  Novelty  Filter  with  Monotony  Reference 

The  value  of  5  (novelty  filter)  was  calculated  from  I(y)/I(x)  =  (1  -  cos5)/(p  +  qcosS). 
p  and  q  needed  to  be  measured  only  once  (see  Appendix  A).  We  measured  I(y)/I(x)  for  100  laser 
shots  and  averaged  them  to  calculate  the  final  value  of  5.  6  can  also  be  calculated  from  I(y)/I(x) 
on  single-shot  basis.  a(5)  is  the  standard  deviation  of  the  single-shot  5  values  in  these  100  shots 
from  the  final  6,  as  given  below. 
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Data  for  Fast  Phase  Modulation  (Novelty  Filter  with  Monotony  Reference) 

In  our  measurement  of  5  (novelty  filter),  we  scanned  the  phase  modulation  by  shifting  its 
time  (relative  to  the  laser)  stepwise  with  1-ns  step  sizes.  At  each  step,  we  took  10  laser  shots  and 
found  I(y)/I(x)  through  a  linear  fit.  After  the  whole  scan,  we  got  a  curve  of  I(y)/I(x).  This 
procedure  was  repeated  10  times  and  the  resulting  10  curves  were  averaged  to  get  a  final  curve. 
From  this  curve  of  I(y)/I(x)  we  calculated  a  curve  of  6  (novelty  filter),  as  shown  in  Figure  6. 

The  phase  modulation  pulse  was  scanned  in  a  similar  manner  during  calibration.  In  the 
measurement  of  5  (novelty  filter)  and  5  (calibration),  the  fast  phase  modulation  was  applied  to 
the  echo  and  the  first  laser  pulse,  respectively.  To  compare  the  two  measured  curves  of  5,  we 
estimated  the  time  interval  between  the  echo  and  the  first  laser  pulse.  We  shifted  the  time  axis  of 
one  curve  of  5  with  respect  to  the  other  and  got  a  reasonably  good  fit,  as  shown  in  Figure  6. 

Discussion 

From  the  data,  we  consistently  saw  a  background  signal  that  gave  a  false  measurement  of 
5  of  15°-20‘’.  This  was  the  result  of  the  incomplete  cancellation  of  the  two  echoes  at  the  output 
when  no  phase  modulation  was  applied.  This  incomplete  cancellation  could  happen  because  the 
two  echoes  at  the  output  were  not  exactly  equal  in  intensity.  However,  the  two  echoes  did  not 
differ  more  than  10%  in  intensity  in  our  experiment.  This  small  intensity  difference  would  not 
give  a  background  in  5  more  than  3°.  Therefore,  the  main  cause  of  background  is  that  the  two 

echoes  are  not  180°  out  of  phase. 
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One  explanation  was  that  the  beam  divergences  of  k2.3  were  not  zero  due  to  diffraction. 
Since  k2  +  k3  had  some  deviation  around  zero,  ke  also  had  some  deviation  around  -ki.  This 
affected  the  phase-conjugate  reflection  from  canceling  the  phase  of  input  laser.  However, 
because  k2,3  are  1-2  mm  in  diameter  and  ki  is  0.2  mm,  we  found  that  diffraction-limited  beam 
divergence  would  not  give  a  background  of  5  more  than  2°. 

The  background  could  also  be  affected  by  how  flat  the  wavefronts  of  k2,3  were  as  well  as 
by  how  flat  the  rare  earth  crystal  surface  was.  We  found  that  the  background  was  quite  sensitive 
to  these  conditions.  The  background  of  5  of  15°-20°  seen  in  our  experiment  could  be  suppressed 
to  about  10°  by  polishing  the  crystal  surface  to  higher  degree.  The  remaining  cause  of  the 
background  needs  further  investigation. 

TEST  OF  PERFORMANCE  ON  MICROSECOND  TIME  SCALE 
Single-Pixel  Phase  Modulator 

We  scanned  the  time  of  the  third  laser  pulse  (relative  to  the  phase  modulation  pulse) 
stepwise,  with  one  step  for  each  single  laser  shot.  The  data  are  shown  in  Figure  7  together  with 
the  calibration  data.  We  shifted  the  time  axis  of  one  curve  of  6  with  respect  to  the  other  and  got  a 
reasonably  good  fit,  as  shown  in  Figure  7. 

In  Figure  7,  the  data  preceding  the  rising  edge  was  the  background,  which  is  about  0.025 
wavelength  (-0.6  |j.m),  or  9°.  The  accuracy  of  the  novelty  filter  was  0.02  wavelength  and  the 
resolution  was  also  about  0.02  wavelength,  or  7°. 

Real  Turbulent  Jet 

The  data  are  shown  in  Figures  8  and  9.  We  used  Eu  doped  crystal  for  these 
measurements.  First,  we  turned  off  the  jet  and  measured  the  background  over  a  long  response 
time  (see  the  lower  parts  of  Figures  8  and  9).  When  the  response  time  was  increased  from  20  to 
200  |j,s,  the  background  clearly  rose.  It  was  0.04  wavelength  at  100  jxs  and  0.08  wavelength  at 
200  )J,s.  We  attribute  this  rise  of  background  to  slow  phase  modulations  along  the  optical  path 
caused  by  the  environment.  These  slow  perturbations  were  rejected  by  the  novelty  filter  under 
short  response  times.  When  the  response  time  was  lengthened,  the  filter  relaxed  to  lower  cutoff 
frequency.  As  a  result,  more  slow  perturbations  leaked  through  the  filter. 
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When  the  helium  jet  was  turned  on,  the  phase  distortion  was  a  continuous  and  random 
event  in  time.  However,  the  novelty  filter  was  limited  to  10-Hz  sampling  frequency  by  the 
laser’s  repetition  rate.  Each  data  point  in  Figures  8  and  9  was  taken  from  the  average  of  100 
sampling  measurements.  The  response  time  of  the  filter  was  then  changed  to  go  to  the  next 
point. 

When  the  laser  crossed  the  jet  at  4  mm  above  the  nozzle,  an  unexpected  sinusoidal 
oscillation  appeared  in  the  novelty  output  as  we  scanned  the  filter’s  response  time  (see  Figure  8). 
The  oscillation  frequency  was  about  10  kHz  and  increased  at  higher  jet  speed.  We  believe  this 
oscillation  is  caused  by  the  jet's  acoustic  vibration,  which  caused  the  jet  thickness  to  change 
periodically  and  induced  oscillation  in  the  optical  path  length.  Although  the  10-kHz  vibration 
was  sampled  at  10  Hz  by  the  laser,  the  periodic  oscillation  can  be  observed  by  averaging  the 
novelty  output.  This  average  is  derived  in  Appendix  B. 

The  jet's  acoustic  vibration  was  proved  by  detecting  the  jet’s  whistle  and  Fourier 
transforming  the  audio  signal.  We  got  the  same  oscillation  frequency  shown  in  the  novelty  filter 
output.  This  acoustic  vibration  is  an  “edge  tone”  of  the  jet  [33].  Its  frequency  can  be  estimated 
by  the  ratio  of  the  jet  speed  to  the  nozzle  diameter.  Using  a  jet  speed  of  100  m/s  and  a  nozzle 
diameter  of  2  mm,  we  get  50  kHz,  which  is  on  the  order  of  the  observed  frequency. 

The  phase  change  was  quite  small  near  the  nozzle  (see  Figure  9)  because  the  jet  in  that 
region  was  laminar  and  changing  slowly  in  time.  The  optical  path  length  was  almost  constant 
there.  Further  away  from  the  nozzle,  the  phase  change  increased  dramatically,  indicating  that  the 
jet  in  that  region  was  turbulent  and  changing  rapidly  in  time.  In  Figures  8  and  9,  the  phase 
changes  rose  to  asymptotic  levels  within  20-40  |xs.  The  rise  time  was  shorter  at  higher  jet  speed; 
40,  30,  and  20  |J,s  at  jet  speeds  of  100, 130,  and  200  m/s.  These  data  indicated  that  the  inverse  of 
the  rise  time  was  linear  to  the  jet  speed.  This  result  is  consistent  with  the  Greenwood  frequency 
[12]  mentioned  above.  As  a  simulation  of  airborne  turbulence,  the  helium  jet  showed  that  the 
phase  distortions  do  change  with  a  bandwidth  of  about  40  kHz. 

In  the  measurements  of  the  jet's  spatial  structures,  we  used  Pr-doped  crystal.  Because  the 
jet  is  a  nonuniform  phase  screen,  the  laser  passing  through  it  and  going  to  the  crystal  directly 
would  be  highly  distorted  at  the  crystal  surface.  As  a  result,  the  laser  would  not  strike  the 
effective  area  of  the  crystal  properly.  This  would  affect  the  phase-conjugate  reflection.  To  solve 
this  problem,  we  imaged  the  jet  on  the  crystal  using  the  laser.  This  image  was  expected  to  be 
uniform  in  intensity,  because  the  jet  was  illuminated  to  be  a  phase-only  object  with  no  amplitude 
modulation. 
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After  we  imaged  the  jet  on  the  crystal,  when  the  novelty  filter  response  time  was  100  ns, 
we  found  that  the  echo  coming  back  through  the  Jet  maintained  a  uniform  profile  at  far  field. 

This  proved  the  effectiveness  of  the  imaging.  The  jet  also  showed  no  phase  change  over  a  shon 
time  of  100  ns.  This  result  indicated  there  was  no  faster  dynamics  in  the  turbulent  jet. 

To  resolve  the  phase  change  on  a  microsecond  time  scale,  we  lengthened  the  response 
time  of  the  novelty  filter  into  a  lO-ps  region.  The  echo  returned  and  illuminated  the  jet  to  be  a 
phase-only  object  that  contained  a  profile  of  the  phase  change  over  the  response  time.  The  jet 
was  then  imaged  on  a  CCD  camera  using  the  echo. 

This  "phase-only"  image  is  shown  in  Figure  10  by  its  intensity  profile.  We  found  that 
the  image  was  not  scattered  and  nearly  uniform  in  intensity.  Inside  the  image  were  some 
nonuniform  structures,  probably  caused  by  some  high  spatial  frequencies  lost  in  the  imaging.  To 
transform  the  phase-only  image  to  an  intensity  profile,  the  echo  carrying  the  image  was 
interfered  with  the  reference  echo.  The  resulting  interference  is  shown  in  Figure  11.  These  data 
showed  the  intensity  of  the  phase  change  over  a  time  of  20  ps. 

We  also  imaged  the  phase  change  over  a  time  of  2  ps  (data  not  shown)  and  observed 
much  weaker  novelty  signal  by  the  interference  (the  turbulence  had  barely  evolved  within  2  ps.) 
We  finally  translated  the  jet  in  x  and  y  directions  transverse  to  the  z  direction  of  the  laser  beam, 
and  we  were  able  to  observe  the  variations  of  the  phase  change  structure  in  different  regions  of 
the  jet. 

ADVANTAGES  AND  LIMITATIONS  OF  THE  PHOTON  ECHO  NOVELTY  FILTER 

The  photon  echo  novelty  filter  has  several  limitations.  So  far,  the  best  sensitivity  of 
photon  echo  novelty  filter  to  phase  change  is  0.025  wavelength.  This  low  sensitivity  is  inferior 
to  the  Shack-Hartmann  wavefront  slope  sensor,  which  has  a  sensitivity  of  0.0006  wavelength 
[14].  The  novelty  filter  has  no  sensitivity  to  the  sign  of  phase  change,  which  means  it  cannot  tell 
a  positive  phase  change  from  negative  change,  because  the  novelty  filter  measures  cos6.  The 
sinusoidal  nature  of  cos6  also  introduces  uncertainty  to  the  absolute  value  of  6  when  5  exceeds 
0.5  wavelength. 

The  photon  echo  novelty  filter  can  be  compared  with  the  traditional  double-pulse 
holographic  interferometer  [26,27].  In  the  latter,  the  phase  distortions  are  recorded  twice  at  two 
consecutive  moments  in  the  same  holographic  film.  After  the  film  is  developed,  the  two 
holograms  are  read  out  simultaneously  and  their  interference  shows  the  time-differential  phase 
change.  Therefore,  the  photon  echo  novelty  filter  and  the  double-pulse  holography  detect  the 
same  change  in  phase. 
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Figure  10. 
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Figure  11. 
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The  key  difference  between  the  two  techniques  is  whether  the  time-differentiation  of 
phase  happens  in  real  time.  The  photon  echo  novelty  filter  needs  to  record  only  one  phase  map 
at  the  first  moment.  This  phase  map  is  immediately  read  out  at  the  second  moment  and  the  time- 
differentiation  is  finished  in  real  time.  This  is  a  great  advantage  of  novelty  filter  over  double¬ 
pulse  holography. 

The  two  techniques  are  different  in  their  data  representations.  The  photon  echo  novelty 
filter  gives  a  map  of  the  intensity  of  phase  change  with  zero  background.  In  contrast,  the  double¬ 
pulse  interferometer  gives  a  map  of  straight  fringes  as  background.  The  amplitude  of  phase 
change  can  be  measured  from  the  distortion  of  these  straight  fringes  [2,3].  Because  it  does  not 
measure  cos5  but  5  directly,  the  double-pulse  holography  has  an  advantage  over  novelty  filter. 

Finally,  the  photon  echo  novelty  filter  is  based  on  inefficient  third-order  nonlinear 
radiations.  The  filter  requires  high  power  lasers  and  cryogenic  cooling  of  the  nonlinear  crystal. 
These  are  drawbacks  compared  with  double-pulse  holography. 


CONCLUSIONS  AND  RECOMMENDATIONS 


In  summary,  we  have  designed  the  setup,  demonstrated  the  concept,  and  tested  the 
perfoiTnance  of  the  photon  echo  novelty  filter  to  measure  a  fast  change  in  phase  distortion.  The 
filter’s  sensitivity  is  better  than  0.05  wavelength  (-0.6  pm)  under  a  response  time  of  <20  ps.  Its 
accuracy  and  resolution  are  both  approximately  0.02  wavelength.  The  filter  is  proved  to  be  able 
to  follow  and  resolve  a  fast  phase  change  event  on  airborne  turbulence  time  scale.  This  novelty 
filter  can  now  be  applied  to  probing  a  table-top  real  turbulent  flow  that  simulates  airborne 
turbulence.  We  showed  that  the  filter  can  resolve  the  coherent  temporal  structures  in  the  flow. 
Spatial  resolution  of  the  filter  was  also  tested. 

In  the  future,  we  wish  to  investigate  several  issues.  First ,  we  have  so  far  proved  that  the 
time-domain-stimulated  photon  echo  can  record  and  recall  a  phase-only  map  through  phase- 
conjugate  recall.  In  other  words,  we  have  proved  that  the  recalled  phase  map  exactly  cancels  the 
original  phase  modulation.  It  would  be  helpful  to  provide  additional  proofs  of  the  phase-only 
memory  by  two  methods.  One  method  is  to  recall  the  conjugate  phase  map,  but  not  let  it  return 
through  the  original  phase  modulator.  Instead,  we  can  use  a  beamsplitter  to  direct  it  to  a 
Hartmann  wavefront  sensor  and  read  it  out.  In  this  way,  we  can  have  a  direct  comparison  of  the 
recalled  conjugate  phase  map  and  the  original  one.  The  second  method  is  to  recall  the  memory 
in  the  forward  echo  rather  than  the  backward  phase-conjugate  echo.  To  do  this,  we  need  to  let 
the  second  laser  pulse  (instead  of  the  first  pulse)  in  the  echo  experiment  carry  the  phase 
modulation.  We  also  need  to  let  the  first  and  third  laser  pulses  be  colinear,  i.e.,  ki  =  ks.  In  this 
way,  ke  =  k2.  We  can  thus  direct  the  echo  to  a  Hartmann  sensor,  read  it  out,  and  compare  it  with 
the  original  one. 

SRI  International  has  been  leading  the  state-of-art  in  amplitude-only  bitmap  memory 
[24].  Based  on  our  understanding  of  the  basic  physics  and  our  expertise  on  the  technology,  we 
believe  a  forward  echo  experiment  will  give  a  definitive  proof  of  the  phase-only  memory  in  hole- 
burning  materials.  We  would  like  to  conduct  such  an  experiment  using  the  turbulent  helium  jet. 
By  directing  a  copy  of  the  laser  carrying  the  original  random  phase  modulation  and  the  echo 
carrying  the  recalled  phase  into  a  Hartmann  sensor,  we  will  be  able  to  have  a  direct  measurement 
of  error  rate. 

We  also  propose  to  apply  the  photon  echo  novelty  filter  to  tracking  the  temporal  structure 
of  fast  phase  change  in  real  time.  To  do  this,  we  need  to  sample  the  phase  change  rapidly 
throughout  the  entire  coherent  time  of  turbulence.  We  would  like  to  acquire  a  laser  source  with 
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high  repetition  rate  (preferably  100  kHz),  although  the  train  of  laser  pulses  need  not  last  longer 
than  the  coherent  time  of  turbulence.  The  laser  should  be  able  to  tune  on  resonance  with  certain 
transition  of  rare  earth  ions. 

We  need  to  modify  the  novelty  filter  so  that  it  can  discern  positive  phase  changes  from 
negative  ones.  One  solution  is  to  install  a  fast  beam  deflector  in  the  reference  arm  of  the  novelty 
filter.  This  beam  deflector  should  be  able  to  provide  a  sudden  and  small  transverse  shift  to  ke 
along  the  reference  arm.  As  a  result,  the  two  echoes  interfering  at  the  output  will  show  straight 
fringes  as  background.  This  will  allow  a  direct  measurement  of  6. 

To  put  the  novelty  filter  on  a  plane  and  probe  real  airborne  turbulence,  the  laser  beam 
diameter  in  the  filter  needs  to  be  expanded  to  about  100  mm  to  suit  the  turbulence  s  spatial  scale 
[2,3].  Therefore,  a  telescope  setup  to  expand  the  beam  and  a  turbulent  jet  with  100  mm  cross 
section  are  needed  in  future  table  top  experiments. 

Finally,  we  want  to  suppress  the  background  of  the  filter  to  a  much  lower  level.  We  have 
already  proved  that  flattening  the  second  and  third  laser  pulses'  wavefront  and  polishing  the  rare 
earth  crystal  surface  give  a  better  phase-conjugate-mirror.  Reducing  the  environmental 
perturbations  also  helps  to  reduce  the  background.  We  want  to  find  out  how  much  background  is 
related  to  the  intensity  of  the  first  laser  pulse  and  its  alignment  of  incidence  on  the  rare  earth 
crystal.  This  knowledge  may  help  reduce  the  background  even  more. 
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APPENDIX  A 


ANALYSIS  FOR  NOVELTY  AND  MONOTONY  CHANNEL  OUTPUT 


If  the  beamsplitter  in  the  interferometer  is  exactly  50/50,  the  ratio  of  novelty  and 
monotony  I(y)/I(x)  is  proportional  to  (1  -  cos5)/(l  +  cos5).  In  reality,  the  beamsplitter  is  not 
exactly  50/50  and  more  careful  analysis  is  necessary  as  follows. 

The  amplitudes  of  the  echoes  on  the  two  arms  are  E  and  E',  respectively  (see  Figure  4). 
The  beamsplitter  in  the  interferometer  has  reflection  and  transmission  coefficients  of  r  and  t  in 
amplitude,  respectively,  r^  and  t"  are  close  to  50%.  The  amplitude  of  the  echoes  that  finally  fall 
on  the  two  detectors  is  Ey  =  (Er  -  E't)a  for  novelty  output  (channel  y)  and  Ex  =  (Et  +  E'r)b  for 
monotony  output  (channel  x),  where  a  and  b  are  constants  that  include  all  loss  factors.  We 
assume  the  intensities  of  the  two  echoes  are  equal  at  the  novelty  output.  In  other  words, 

I(y,e)  =  I(y,e')  =  1(e)  when  either  arm  is  blocked.  Under  this  condition,  the  two  echoes  at  the 
novelty  output  differ  only  in  phase,  i.e.,  E't/Er  =  e^^. 

Let  f  =  (r/t)^  and  q  =  (b/a)^.  Then  the  novelty  output  is 


I(y)  =  EyEy*  =  2(1  -  cos5)I(e) 


and  the  monotony  output  intensity  is 


I(x)  =  ExEx*  =  2(p  +  qcos5)I(e) 


where  p  =  q(f^  +  l)/2f.  Finally,  I(y)/I(x)  =  (1  -  cos5)/(p  +  qcos5). 

To  find  p  and  q,  we  need  to  measure  the  ratio  between  the  two  channels  when  one  of  the 
arms  is  blocked.  In  the  novelty  output,  I(y,e)  =  I(y,e')  =  1(e).  In  the  monotony  output, 

I(x,e)  =  I(e)(q/f)  when  E'  is  blocked,  and  I(x,e')  =  I(e)(qf)  when  E  is  blocked. 

Thus,  if  we  measure  the  ratio  s(e)  =  I(y,e)/I(x,e)  when  E'  is  blocked  and 
s(e')  =  I(y,e')/I(x,e’)  when  E  is  blocked,  then  f^  =  s(e)/s(e')  and  (1/q)^  =  s(e)s(e').  From  this,  we 
find  p  and  q. 
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APPENDIX  B 


DERIVATION  OF  AVERAGE  OF  NOVELTY  OUTPUT  WHEN 
SAMPLING  PERIODIC  OSCILLATION 

The  novelty  filter  output  is  ~[1  -  cos5],  where  5  is  the  change  of  optical  path  length  over 
the  response  time.  When  5  <  1,  the  filter's  output  is  ~  5^. 

For  a  sinusoidally  oscillating  optical  path  length, 


5  ~  sin[a)(ti  +  T)]  -  sin(coti) 


where  co  is  the  oscillation  frequency,  q  is  the  time  when  the  laser  samples  the  medium,  and  T  is 
the  novelty  filter's  response  time.  Thus 


52  ~  sin2[a)(ti  +  T)]  +  sin^fcoti)  -  2sin[0D(ti  +  T)]sin((Oti). 


Because  the  laser  is  not  synchronous  with  the  oscillation,  tj  is  a  random  sampling  time.  If 
we  take  the  average  of  many  laser  shots, 

<62>  ~  1/2  +  1/2  +  <cos[co(2ti  +  T)]  -  cos(cdT)  >  =  1  -  cos(coT). 
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